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NONMONOTONICITY OF PHASE TRANSITIONS IN A LOSS 
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and University of Auckland 

We consider a symmetric tree loss network that supports single- 
link (unicast) and multi-link (multicast) calls to nearest neighbors 
and has capacity C on each link. The network operates a control 
so that the number of multicast calls centered at any node cannot 
exceed Cv and the number of unicast calls at a link cannot exceed Ce, 
where Ce, Cv < C. We show that uniqueness of Gibbs measures on 
the infinite tree is equivalent to the convergence of certain recursions 
of a related map. For the case Cy ~ 1 and Ce = C, we precisely 
characterize the phase transition surface and show that the phase 
transition is always nonmonotone in the arrival rate of the multicast 
calls. This model is an example of a system with hard constraints that 
has weights attached to both the edges and nodes of the network and 
can be viewed as a generalization of the hard core model that arises in 
statistical mechanics and combinatorics. Some of the results obtained 
also hold for more general models than just the loss network. The 
proofs rely on a combination of techniques from probability theory 
and dynamical systems. 

1. Introduction. In Section 1.1 we describe and motivate the basic loss 
network model analyzed in this paper. In Section 1.2 we summarize the main 
results and place our results in the context of related prior work. An outline 
of the paper is provided in Section 1.3. 

1.1. Description of the model and motivation. Many stochastic processes 
that arise in applications take values in 5^, where 5 is a finite state space 
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and V is the vertex set of a finite or countably infinite loopless graph 
G = {V, E) with edge set E. A generic example is the Ising model of fer- 
romagnetism [1], where S = {—1,1} represents the set of possible magnetic 
spins at a vertex site of G, and configurations in that have opposing 
spins on neighboring sites are discouraged (i.e., assigned lower probabil- 
ity). Another well-known example is the hard core model, which was first 
proposed to study the equilibrium behavior of a lattice gas consisting of par- 
ticles with nonnegligible radii that cannot overlap [1, 14, 23]. Here, the state 
s S 5 = {0, 1} of a node represents the occupation number of a particle, and 
particles cannot occupy neighboring vertices on the lattice. In other words, 
the set of feasible configurations is given by 

nQ' = {aeS^ ■.a^ + ay<iyx,yeV such that xyGE}. 

The hard core model is an example of a stochastic process subject to "hard 
constraints," in which certain configurations are forbidden (rather than just 
assigned lower probability, as in the case of the Ising model) . Processes with 
such "hard constraints" arise in fields as diverse as statistical mechanics, 
combinatorics and telecommunications. In particular, the hard core model 
also arises in the study of random independent sets of a graph [5, 12, 20] and 
in the analysis of multicasting on telecommunication loss networks [18, 26]. 
The consideration of different applications has led to the study of various 
generalizations of the hard core model [2, 5, 13, 22, 24, 26, 28, 29]. One natu- 
ral generalization of the hard core model arises from the study of stochastic 
loss networks in telecommunications. 

A general loss network consists of E links (or resources), indexed by e, 
with the eth link having integer capacity Cg. Calls are indexed by a set TZ, 
and type r calls arrive as a Poisson process with rate jr and stay in the 
network for a random duration. Upon arrival, a call of type r requests ca- 
pacity from the eth link for each e, where the set of links {e : A^r > 0} is 
often referred to as the route of the call. If the requested capacity is available 
simultaneously on all links along its route, then the call is accepted; other- 
wise it is rejected and lost. If the call is accepted, it occupies the requested 
capacity for its duration. The stationary distribution for a loss network de- 
pends on the call duration distribution only through its mean [7] , so we allow 
call durations to have a general distribution with finite mean, which, with- 
out loss of generality, can be assumed to be equal to 1. All arrival processes 
and call durations are assumed to be independent of one another. For an 
excellent introduction to loss networks, see the review paper by Kelly [18]. 

In this work, we consider loss networks with the structure of a graph 
G = {V, E) that is a regular tree (though the description in this section 
holds for a general graph G). All links or edges in the network have the 
same capacity C, and the network supports two classes of calls — multicast 
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calls that arrive at a node with a Poisson rate v and require one unit of 
capacity from each of the edges incident to it, and unicast calls that arrive 
at an edge with a Poisson rate A and occupy one unit of capacity on that 
edge. In addition to the capacity constraints that determine whether or not 
a call can be accepted, we also impose a simple control, whereby a call 
is accepted only if, after acceptance, there will be no more than Cy < C 
multicast calls at any node and no more than Ce unicast calls at any 
link in the network. Note that when Cy = Ce = C , the additional admission 
control is redundant and the controlled unicast-multicast model reduces to 
the uncontrolled unicast-multicast model that was introduced in [26]. 

Let n(s) = {nv{s),ne{s),v €V,e £ E}, where n^(s) and ne(s) represent 
the numbers of multicast and unicast calls in progress at time s at node v 
and on edge e, respectively. Then, given any graph G = {V,E), under the 
above admission rule for every s € (0, oo), n(s) takes values in the set ilc of 
feasible configurations given by 

(1.1) CIg = {w € S^y X S^^ : uJu + i^uv + i^v < C Vti, V such that uv G E}, 
where for any positive integer C, 

(1.2) Sc = {0,l,...,C}. 

Given vectors u = (z/q, z^i, . . . , ucv) ^ and A = (Aq, Ai, . . . , \cj € 

jgCfi+i^ we define mj^'*' to be the following probability distribution on CIg- 

(1.3) m^j (uj) = < z-.^ 

[ 0, otherwise, 

where Zq^ is the usual normalizing constant defined by 

Note that the quantity Hi-ey IleG-B At^e is proportional to the probability 
of the configuration uj. We will refer to z7 and A as weight vectors (they are 
sometimes also called activity vectors in the literature). 

For finite G, it is well known [18] that the process n(-) has a unique 
stationary distribution, which is supported on CIg and is of product-form, 
given explictly by rriQ^, with weight vectors 

I/* 

i/j = — for ? = 0,l,...,Cy, 

(1.4) 

A* 

Aj = — for I = 0, 1, . . .,Ce- 

ii 
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Here, we adopt the usual convention that 0! = 1, so that z^o = ^^o = 1- The 
structure of the weight vectors given in (1.4) is a consequence of the assump- 
tion that calls arrive as a Poisson process and are blocked if they cannot be 
served immediately (there is no queueing), with each accepted call allocated 
one unit of capacity. Thus, the downward rate of calls at a node with occu- 
pancy i is exactly i (since call durations have mean 1). On the other hand, 
networks of interest often support processor sharing rather than allocate a 
fixed capacity to each accepted call. In the latter case, an arriving call is 
accepted only if the buffer capacity C is not exceeded at every link along 
which the call requires capacity (so that the set of feasible configurations 
remains the same), but all multicast calls present at a node are served in a 
processor-sharing manner by a processor sitting at the node and, likewise, 
all unicast calls present at a link are served in a processor-sharing manner 
by a processor at the link. Motivated by the processor-sharing discipline, in 
the treatment below, we will also consider product-form distributions of the 
form (1.3) with the weight vectors z? and A defined instead by 

(1.5) fi = i^* for i = 0,l,...,Cy and Aj = A* for z = 0, 1, . . . , Ce. 

Such weights have also been considered for the uncontrolled pure multi- 
cast model on the lattice in [24]. Another purpose of considering alternative 
weights is to gain some insight into the influence of different weight vectors 
on the existence of phase transitions. 

More generally, we will refer to any model possessing given weight vec- 
tors z7, A and controls Ce, Cy as the controlled unicast-multicast model, 
recognizing that, in the special cases where the weights are given by (1.4) 
and (1.5), these have a specific interpretation in the communications context 
as models for loss and processor sharing models, respectively. 

Although the product measure (1.3) has an explicit representation, for 
large graphs, the exact computation of the stationary probabilities is made 
difficult by the combinatorial complexity of calculating the normalizing con- 
stant. Instead, in such cases, one often studies the behavior of the measures 
in the asymptotic limit as the size of the graph increases. The appropri- 
ate measures on unicast-multicast configurations on the infinite limit graph 
G = (y, E) that provide insight into the behavior of product-form measures 
on sufficiently large finite graphs are Gibbs measures. A Gibbs measure is a 
distribution of a countably infinite family of random variables which admits 
some prescribed conditional probabilities [14]. Roughly speaking, the Gibbs 
measures studied in this work are characterized by the property that the dis- 
tribution of the configuration on any finite subset U oiV\JE, conditioned 
on the configuration on the complement, is equal to the regular conditional 
probability of an associated product-form measure of the configuration on 
[/, given the configuration on the boundary of U (see Definition 3.2 for a 
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more precise formulation). Unlike stationary distributions on finite graphs, 
the associated Gibbs measures on infinite graphs need not be unique. In par- 
ticular, there may be multiple Gibbs measures associated with a given pair 
of weight vectors. A phase transition is said to occur if one can identify the 
boundaries between weights for which there is a unique Gibbs measure and 
weights for which there are multiple Gibbs measures. Moreover, the phase 
transition is said to be monotone with respect to a certain parameter if the 
existence of multiple Gibbs measures at a certain value s of the parameter 
implies the existence of multiple Gibbs measures for all parameter values 
greater than s. (Note that, since one can always consider the reciprocal of 
the parameter, there is no loss of generality in assuming that multiple Gibbs 
measures persist only with an increase in the parameter.) An illuminating 
discussion of the implications of phase transitions on the infinite graph for 
the nature of the stationary distribution on large finite graphs can be found 
in [19]. 

An important motivation for studying this model is to gain a better 
understanding of the general structure of Gibbs measures associated with 
product-form distributions on graphs that have weights on edges in addi- 
tion to nodes. Such measures have received much less attention than those 
where weights are just on the nodes, which have been studied, for example, 
in [5, 13, 17, 28, 29]. The uncontrolled model in the presence of unicast and 
multicast calls (with Cy = Ce = C) for the simple case C = 1 was studied 
in [26]. As discussed in further detail in Section 6, we expect that the con- 
trolled model here (with Cy = 1) may provide insight into the behavior of 
the uncontrolled unicast-multicast model when C > 1. 

Interest in this model initially arose from the need to gain insight into 
the effects of multicasting on the performance of a large network. Multi- 
casting is a major recent innovation in communication networks, which was 
introduced to deal with the emergence of new multimedia applications that 
require simultaneous connections between multiple users (as in audio or 
video conferencing) and simultaneous transmissions to multiple users (as in 
the Internet) [8, 25]. In [22, 26] we showed that models with multicasting 
can give rise to phase transitions. The occurence of such phase transitions is 
undesirable, because customers experience this as unfairness in the network 
(in the form of spatially heterogeneous blocking probabilities, even though 
arrival rates are homogeneous). It is natural to consider whether controls 
can be used to mitigate this effect, and the control we propose in this paper 
is a particularly simple one. It has a similar flavor to traditional trunk reser- 
vation, but because reversibility is preserved, the stationary distribution still 
has product-form (unlike the case of trunk reservation). 

1.2. Main results and related prior work. The first main result of the 
paper. Theorem 1.1, shows that, given a pair of weight vectors A and P, 
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uniqueness of Gibbs measures for the controlled unicast-multicast model 
described above is equivalent to the convergence of certain recursions of 
the associated map ^'^'■^ introduced below. Given positive integers Cy, Ce 



for k = 1,. . . , Cy and G . We refer to the mapping as the ran- 
dom field map associated with the weight vectors z7 and A. We now state 
Theorem 1.1. For a precise definition of Gibbs measures associated with the 
controlled model, refer to Definition 3.2. In the following, denotes the zero 
vector: 

Theorem 1.1. Let ^'^'■'^ be the random field map defined in (1.6). Then, 
the unicast-multicast model on the (q -\- 1) -regular tree, G, has a unique Gibbs 
measure if and only if, given the initial condition ^^^^^ = G M^^, the se- 
quence defined iteratively by = $^^'^(^("-1)) converges to a limit. In this 
case, the limit will be the unique fixed point of the map (^^'^ . 

The proof of Theorem 1.1 is given at the end of Section 3.2. Similar corre- 
spondences have been obtained for the case of Gibbs measures on a feasible 
subset of the configuration space for some finite set 5", when the weight 
vectors are attached to just the nodes f G V of a regular tree (see [27] for 
the case of a binary state space, [28, 29] for the general case of bounded, 
positive Markov specifications on a countable state space and [5] for more 
recent work on the case of product-form Markov specifications). However, 
when the edges also have weights and the set of feasible configurations is 
given by (1.1), then the graph with respect to which the specifications are 
Markov is no longer a tree and, thus, the results of [5, 28, 29] are no longer 
directly applicable. Theorem 1.1 extends this correspondence to the more 
complicated setting in the presence of unicast calls — where the edges also 
have weights, and the specifications are Markov with respect to the aug- 
mented graph G defined in Section 3.1.1, rather than with respect to the 
regular tree G. 

Analysis of recursions for the uncontrolled unicast-multicast model was 
carried out in [26]. In particular. Theorem 1.1 also shows that the results 
of [26] have implications for the corresponding Gibbs measures. Specifically, 
it rigorously proves that when C = 1, the value obtained in (3.7) of [26] is 
indeed the phase transition point for the (uncontrolled) unicast-multicast 




, where 
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model and, also, that phase transitions continue to occur for the hard core 
model on the tree, even when the Markov specifications are nonhomogeneous 
and parity-dependent (see Theorem 5.1 of [26]). The latter property is in 
contrast to the behavior of the hard core model on the lattice, where phase 
transitions disappear when the Markov specifications are nonhomogeneous 
[2, 15]. 

The second important result of this paper is a precise characterization of 
the phase transition surface (namely the boundary between where there is 
uniqueness of Gibbs measures and where there are multiple Gibbs measures) 
for the controlled unicast-multicast model on the q + 1-regular tree with 
Cv = 1, Ce = C and weight vectors A that satisfy Assumption 1.1 stated 
below. Note that, when Cy = 1, since 1^0 = 1, the weight vector z7 can be 
identified with the scalar u = I'l. 

Assumption 1.1. For every C > 2, the weight vector A = (Aq, Ai , . . . , Ac) 
satisfies 

(1.7) Al_^ - AcAc-2 > 0, 

where Ac = Ei=o-^i- 

As shown in Lemma 4.1, this assumption is satisfied, in particular, by 
the weight vectors of interest described in (1.4) and (1.5). The set A, the 
mapping J and the quadratic Q defined below will be useful in characterizing 
the phase transition surface. 

Definition 1.1. We define A to be the subset of the {q, C, A)-parameter 
space {2,3,...}^ x (UdLil^+) for (g, C, A) that is characterized by the con- 
dition that {q, C,X) ^ A if and only if 

{q + lfAcAc-2 <{q- 1)^ Ac-i • (Condition A) 
The mapping J : [0, oo) — s- [0, oo) is defined by 

and the quadratic polynomial Q is equal to 
Q{a) = a^kc-iKc-2 

(1.9) 

+ a[(l - g)A^_i + (1 + g)AcAc-2] + AcAc-i. 

The result can now be summarized in the following theorem, the proof of 
which is given at the end of Section 4.1: 
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Theorem 1.2. Given A satisfying Assumption 1.1, the controlled unicast- 
multicast model with Cy = 1 and Ce = C on the {q + 1) -regular tree associ- 
ated with the parameters {q, C, A, z^) has multiple Gihhs measures if and only 
if (?) C,X) £ A and v £ z^+), where u- = J{a-) and iy+ = J{a+), with J 
being the mapping defined in (1.8) and a_ and a+ being the distinct positive 
real roots of the quadratic Q specified in (1.9). 

Note that the above theorem shows that whenever the phase transition 
occurs, it is nonmonotone in the parameter u, for fixed (g,C, A). In general, 
determining monotonicity of Gibbs measures is a nontrivial task. Indeed, the 
question of monotonicity of phase transitions for even the hard core model 
on certain graphs such as the cubic lattice remains an open problem. 
Interesting recent results have shown that the hard core model is monotone 
on certain n-dimensional lattices [16], but nonmonotone on certain other 
graphs [4]. The graph on which nonmonotonicity was demonstrated in [4] is 
a somewhat artificially modified tree, having several pendant nodes hanging 
from each node. In our work, nonmonotonicity is shown to occur in a related 
model arising in the context of a real application. 

1.3. Outline of the paper. The outline of the paper is as follows. In Sec- 
tion 2 we relate the random field map 

to recursions that characterize 
the limiting stationary distribution at a central node of a large finite tree 
network and obtain expressions for the blocking probabilities. In Section 3 
we provide some general background on Gibbs measures, establish proper- 
ties of the Gibbs measures associated with the controlled unicast-multicast 
model and present the proof of Theorem 1.1, which reduces the study of 
uniqueness of Gibbs measures to the analysis of recursions of the random 
field map. In Section 4 we analyze these recursions and study the phase tran- 
sitions for the case Cy = 1 and Ce = C. Numerical results are presented in 
Section 5 and a concluding discussion is given in Section 6. 

2. Limiting stationary distributions. In this and subsequent sections, we 
only consider the case when the graph G has the form of a {q -\- l)-regular 
tree T [recall that a {q -\- l)-regular tree is the unique graph with no cycles 
such that each node has precisely q -\- 1 neighbors]. Although the compu- 
tation of the normalizing constant Zq^ and, therefore, of the stationary 
distribution, is hard for general large graphs, the calculation is somewhat 
simplified when the graph is a tree. In particular, due to the absence of cy- 
cles in a tree, it is possible to derive a recursion relation for the normalizing 
constants for trees of increasing size (see [26] for the recursion relation for 
the uncontrolled version of this model) . 
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Let Tm be the complete g-ary tree of height m, with the root node denoted 
by O. In other words, is a rooted tree in which the root node, O, has 
degree q, the qm terminal or leaf nodes have degree 1 and all other nodes 
have degree g + 1. Fix z? and A and let Qm = ^T„i- In order to derive the 
recursion relation, we take advantage of the property that the removal of any 
edge in a tree splits that tree into two disjoint trees. For i G {0, 1, . . . , Cy}, 
let Zm{i) be the weighted sum of all feasible configurations on C T that 
have i multicast calls at the root node O. More precisely, define 

Zm{i)= Yl^n^Il^n, 

Then, the vector {Zm.^i{i),i = 0,1, . . . ,Cv) for the tree Tm+i can be ex- 
pressed in terms of that for as follows: 

(min(C— ijCg) min(C—i—j,Cv) \ 1 

j=0 k=0 ) 

Since a finite (g + l)-regular spherical tree T (say of diameter 2L) can be 
decomposed into a graph containing the central node with g + 1 edges 
incident to it and g + 1 rooted trees each of height L — 1, the normalizing 
constant for the tree T is given by 



i=0 



where 



(min(C— ijCs) min{C-i— j,Cy ) \ 
Y Y ZL^i{k)\ 

j=0 k=0 ) 

Given boundary conditions for the external nodes, the normalizing con- 
stant Z)^^ can be calculated from the recursion (2.1) and the expressions in 
the last two displays. Thus, the problem of computing the stationary distri- 
bution of a finite tree is reduced to that of analyzing the recursion relation 
(2.1). In Section 3 (see the proof of Theorem 1.1) we show that the recursion 
relations are also useful in determining the structure of Gibbs measures on 
the infinite (g -|- l)-regular tree. 

The components of the vector {Zm{i),i G {0, 1, . . . ,Cv}) tend to infinity 
as m — > oo. However, the stationary probability that the occupancy at the 
root of Tm is equal to i, where i € {0, 1, . . . , Cy}, is bounded and lies in 
[0,1]. Thus, it is more convenient to consider a renormalized vector, which 
carries adequate information for computing the stationary distribution. For 
i = 0, . . . , Cv and m e N, define 

^ Zmji) 
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and let = (^^(1), • • • ,(,m{Cv))- Note that is well defined since ^m(O) 
is strictly positive for every m. The recursion (2.1) can then be recast into 
the more convenient form 

where, as defined in (1.6), is the random field map associated with the 
weight vectors u and A. 

For a loss network, important performance measures of interest are the 
blocking probabilities of different types of calls. Since calls arrive as Poisson 
processes, the stationary blocking probability of a call is equal to the sta- 
tionary probability that the system is in a configuration state such that the 
addition of the call to the network would take the configuration out of the 
feasible set CIq. For instance, if (3^ is the blocking probability of a multicast 
call arriving to node v on the tree Tl, we can write 

Y^C'v--! |'Y^min(C-i-l,Ci5) , ^min{C-i-j-l,Cv) c ^ ,^^q+l 

An expression for the blocking probabilities of unicast calls can also be de- 
rived in a similar fashion. From these expressions, it is possible to deduce 
that, in contrast to multicast blocking probabilities, the unicast blocking 
probabilities are always spatially homogeneous, in the sense that they de- 
pend only on the distance from the center of the network and, in the limit 
as L — > cxD, the links at the center of the network have the same blocking 
probability. 

3. Gibbs measures. In this section we introduce Gibbs measures for the 
controlled unicast-multicast model on the infinite {q + l)-regular tree de- 
scribed in Section 1.1. Section 3.1.1 introduces some basic notation and 
Section 3.1.2 provides general background on Gibbs measures and Markov 
specifications. Section 3.2 focuses on Gibbs measures associated with the 
controlled unicast-multicast model, with the main results being Theorem 
3.1 and the proof of Theorem 1.1. 

3.1. Notation and definitions. 

3.1.1. Basic notation. In this section, let G = {V,E) be the graph rep- 
resenting a regular infinite tree network with nodes or vertices V and links 
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or edges E. We will occasionally write f ~ to if vw G E. Let G = {V,E) be 
the augmented graph defined by 

V = VUE, 

E = EL) {xe : x E y, e G -E and xz = e for some z € V}. 

For any U CV ,\et dU = {x eV\U : xz e E for some zeU} and U = UUdU 
denote, respectively, the boundary and closure of U with respect to the 
neighbor relation in E. Likewise, for any U CV, let dU = {x (^V \ U : xz 
E for some z £ U} and AU = U U dU denote, respectively, the boundary 
and closure of U associated with the neighbor relation in E. For U cV, for 
conciseness we will often use Uy to denote U CiV and Ue to denote U (1 E. 
We will also write V (resp., E) for the set of nodes in V that correspond 
to nodes (resp., edges) in G. In general, ii AcV, we will let ^ C F denote 
the set of nodes in V corresponding to the nodes in ^4 C V^, and likewise for 
B C E. For U cV, we use G[U] to denote the graph with the set of vertices 
equal to U and the set of edges comprising xy £ E, such that x,y £U. For 
U cV, G[U] is defined analogously. 

For U C V , au : S^^ — > is the projection <Jij{u}) = {iOi,i € U} onto U, 

and cry is simply denoted by a. Analogously for U cV, au -Sq^ x Sq^ — > 

X is the projection au{uj) = {u!i,i £ U} onto U and, again, the 
identity mapping ay is denoted simply by a. Moreover, for convenience, cjj 
and ai will be used instead of ajj} and a^ij , respectively. Let V and V be the 
set of finite nonempty subsets of V and V, respectively. For U,W cV with 
UnW = 0, T£S^^ andr]£S^^, tt] = {t^ : x £ U} U {rj^r : x £ W} £ S^^^^ 
represents the concatenation of r and r], and the obvious analogue holds 
for configurations in Sq^ x S^^. For U cV, let ^{U) be the u-field in Sq^ 
generated by sets of the form {ai = r} for some i £U and r £ Scy ■ Likewise, 
for UCV, let be the a-field in S^"^ x 5g generated by sets of the 

form {a-i = r}, where either i £ Uy and r £ Scy , or i £ Ue and r £ Sqe ■ 

For conciseness, henceforth in this section we will use S to denote Scy 
and, for U cV, we will use 5^ to denote S^^ x S^'^. The set given in (1.1) 
of feasible configurations for the controlled unicast-multicast model on any 
graph G with vertex set V and edge set E can then be rewritten as 

(3.1) CIg = {oj £ : {j)x + '^xy +ujy <G for all x,y £V such that xy £ E}. 
Note that if we define 

(3.2) 0,G = {uj £ -.ijJx + u}y <G for all x,y £V such that xy £ E}, 
then, for every lo £ (Ig, o"y(cj) £ Q.g- 
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3.1.2. Markov specifications and Gibhs measures. We recall here the defi- 
nition of a Markov specification (see, e.g., [14, 28]), adapted, where necessary, 
to the context of this application. 

Definition 3.1. A Markov specification on {S^ ,T{V)) is a collection 
n = {■^f/}(7gv of stochastic kernels : S'^^ x T{U) —>■ [0, 1] that satisfy the 
following consistency condition: for each U,W £V with U CW and r S , 

. , . . T^ui^Quiau = Tu) 

^w{rQyy;(Tw = Tw) = ^ ; ^ 

'I'W,u[taw\u) 

for some normalizing function hw,u '■ S^^'^^ — > M+. 

Recall that, by the definition of a stochastic kernel, for each U £V and 
T € , ttu{tqjj] •) is a probability measure on {S^ ,!F{y)). 

A generic way of obtaining Markov specifications on an arbitrary graph 
G is to consider the collection of regular conditional probabilities associated 
with a consistent family of "product-form" probability distributions defined 
on all finite subgraphs of G. Specifically, given (weight) vectors v € R'-'^"'"^ 

and A G R^^"*"^, for any [/ E V, first let the associated distribution rhu = 
be defined on by 

(3.3) rriu{a = T) = \ ' , liT^^GiU], 

y 0, otherwise, 

where Zjj is the usual normalizing constant that makes rhu a probability 
measure and ^g[U] is defined by (3.1). Note that this measure is precisely the 
stationary distribution on the finite graph G\U] of the controlled unicast- 
multicast model associated with weight vectors A and z?, as described in 
Section 1.1 and, in particular, for the weights given in (1.4). Then, for [/ E V, 
let TTf/ : S^^ X T{U) be the stochastic kernel defined by 

(3.4) t^u{tqu',^w = tw) =rhAu{^w = '^w\^qu = Tdu) 

for every r^^ G with m/:^u{crQjj = Tqjj) > 0, W QU and tw € It IS 
straightforward to check that this defines a consistent family of stochastic 
kernels in the sense of Definition 3.1. 

While the expression in (3.3) is well defined for finite sets U, it does 
not make sense for infinite sets, since Zjj would be infinite. The natural 
generalization of the stationary distributions {?wc/}f/gy to infinite graphs 
turns out to be Gibbs measures that correspond to the Markov specification 
n = {vrc/}c/ev, with % given by (3.4). 
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Definition 3.2. Given a Markov specification tl on (S"^, JF(V)), a Gibbs 
measure corresponding to IT is a probability measure jl on {S^ ,T{V)) that 
satisfies, for all C/ G V and fi a.s. r € S"^, 

(3.5) fi{au = Tu\&y\jj = Tyy^) =Tfu{Tg^;au = Tu). 

A Gibbs measure fi that corresponds to the particular Markov specification 
n = {'^u}ij^y defined by (3.4) will be called a Gibbs measure for the con- 
trolled (unicast-multicast) model. We let G{Il) denote the space of Gibbs 
measures corresponding to the Markov specification 11 . 

In Section 3.2 (see Corollary 3.7) we establish a one-to-one correspondence 
between Gibbs measures for the controlled model and certain Gibbs mea- 
sures on (5^,J^(F)). For this, we need to introduce analogous definitions of 
Markov specifications and Gibbs measures on . 

Definition 3.3. A Markov specification on {S'^ ,J^{V)) is a collection 
n = {t^u}u&v of stochastic kernels ttu : S^^ x !F{U) [0, 1] that satisfy the 
following consistency condition: for each U,W £V with U C W, and r G S^: 

, N '^u{Tdu;cru =Tu) 

T^w[Taw](yw = Tw) = — r t ^ 

n'W,u[TAw\u) 

for some normalizing function hw,u '■ S^^'^^ — > M+. 

Definition 3.4. Given a Markov specification U = {nuju^y on (5^, 
J-{V)), a Gibbs measure corresponding to IT is a probability measure fi on 
{S'^ ,J-{V)) that satisfies, for all [/ G V and /i a.s. r G , 

(3.6) fi{au = Tu\crv\u = tv\u) = t^u{tqu;cfu = tu)- 

Analagously to Definition 3.2, we will denote by Qiji) the space of Gibbs 
measures corresponding to the Markov specification 11 . 

It turns out that the Markov specifications on {S^ ,J-iy)) that arise in 
our analysis cannot be generated from a family of consistent product-form 
probability distributions. Instead, we will require a more general way of 
generating a Markov specification on {S^ ,J-{y)), namely, via an interaction 
function. Specifically, a mapping (p:SxS^ M+ is said to be an interaction 
function if and only if there exists some reference element /cq G 5 such that 

(/)(i, ko) > 0, i G 5, 

<t>{iJ) = 4>{3,i), hj^S. 

(We have omitted a third condition on the mapping cf) given in [28] that 
is needed to ensure that the associated Markov specification is bounded. 
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since it is automatically satisfied here, due to the finiteness of S.) Given an 
interaction function (j), it was shown in [10, 11] that the collection {'n'(f,^u}ueV 
given by 

defines a consistent family of stochastic kernels tt^^jj : S'^^ x J-{U) M_|_, 
in the sense of Definition 3.3. Here k^^jj : —)■ R+ is the appropriate nor- 
malizing mapping that, for each rgu € S^^ , makes TTip^ui^dU] ") ^ probability 
measure. Given an interaction function (p and the collection ir^^jj, [/ S V, de- 
fined by (3.7), we use 11^ = {TT(j,^u}ueV to represent the Markov specification 
with interaction function (p. 

3.2. Gibbs measures for the controlled model. The first main result of 
this section is Theorem 3.1 which, for a given pair of weight vectors z7 and 
A, establishes a one-to-one correspondence between the set G{tl) of Gibbs 
measures for the controlled model and the set ^(H^) of Gibbs measures 
corresponding to the Markov specification with interaction function cp 
defined by 

mm(CE,C-{i+j)) 

(3.8) Hhj) = { k=o 

for i,j e S:0<i+j <C, 

0, otherwise. 

Note that this function (p is clearly symmetric and serves as a reference 
element, so that (p is indeed an interaction function. In fact, it gives rise to 
a repulsive Markov specification, in the sense defined by Zachary [29], since, 
for all i <:k,j <l, 

4>{hj) < 4>{k,j) 



<p{i,l) - (p{k,l) 

The proof of the theorem relies on the following two lemmas. The first lemma 
is an elementary graph-theoretic result, while the second lemma establishes 
a relation between the specifications 11 and II,^. 

Lemma 3.5. For any U CV, d[U U [£■ n dU]] = dU . Moreover, for any 
U dV, d{AU nE) = Aunv. 

Proof. First, note that for U CV CV, by the definition of d and d, it 
follows that dU=[En dU] U dU and d[E ndU] C dU UU . The first relation 
in the lemma then follows from the observation that 

d[A UB] = [dA U dB] \[AuB] for all A,BcV. 
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For the proof of the second relation, first note that, by the last display 
and the fact that Ac E implies dA C V, we have d[AU D E] = dpE U 
[dUv n E]] = BUe U d[dUv n ^] and AC/ n y = BUe U Uv From the last 
two equalities, it is clear that to complete the proof, it suffices to show that 
Uv = d[dUy n E]. Now, suppose x ^Uy. Then, there must exist 2: G C/y 
such that xz G dUy H E, from which it directly follows that x G d\dUv H E]. 
On the other hand, if x G d[dUv n E], then there must exist y G dUy n E 
such that y = xz for some z €V. However, the fact that y G dUy H E implies 
that y = uv for some u G Uy and v G C/y . Since the representation for y is 
unique, this implies that x = u or x = v, from which it follows that x (^Uy, 
and the proof is complete. □ 

Lemma 3.6. Let H be the Markov specification for the controlled unicast- 
multicast model and let II^ be the Markov specification associated with the 
interaction function 4> defined by (3.8). Then, for any r G and U CV , 

'^Uu[EndU]('^9u;^U = Tu) =T^<t,,u{rdu;cru = Tu)- 

Proof. Fix r G and let U cV. Then, by Lemma 3.5, dU = d[U U 

[E n dU]], and so TT^u[EndU](^du;^u = w) = T^uu[EndU](^d[Uu[EndU]]'^^u = 
m) which, by (3.4) and the fact that A[UU[En dU]] = UiJ[EndU]UdU = 
AC/, is equal to 

iriAui^u = ^ul^du = Tdu) 

= ^^ui^uu[dunE] = -ruvl^du = Tdu) 

rhAui^dU = Tdu) 

nTT TT •^rain{CE,C—Tx—Ty) , 
_ u(^U liv&dU llxyGAUnE l^k=0 ^^i 

mAui^du = Tdu) 

if < + < C for uv (z E f] AC/, and otherwise. For v G dU, let = 
\{uv:u (zU,v dU}\ be the number of neighbors of f in C/ (clearly r^ < 
q + 1). Then, each u^U has g + 1 edges in AC/ emanating from it and 
each V G dU has r„ edges in AC/ emanating from it. Consequently, when 
< Tu + Tt, < C for every uv G Ef] AC/, the right-hand side of the last display 
can be rewritten as 

ruAui^dU =Tdu) 
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Now, define 

\vedu ) 

Then, combining the last three displays and recalling the definitions of 
and Tx^^u given in (3.8) and (3.7), respectively, we obtain 

^uuiEndU]'^^du;(^u = Tu) = k^u{TQu) " ^^^uijau-^'yu = Tu), 

which completes the proof of the lemma. □ 

Given a probability measure fi on , the projection n oi jl onto V is 
defined in the obvious way: for all r G and U QV, 

fi{au = Tu) = fi{uj G -.auiuj) = m) = K^u = m)- 
We now state and prove Theorem 3.1. 

Theorem 3.1. Suppose positive integers Cv,Ce,C and weight vectors 
are given. Let 11 he the Markov specification for the 
associated controlled unicast-multicast model and let 11^ he the Markov spec- 
ification with interaction function (j) defined hy (3.8). Then, the following 
three properties hold: 

1. Given fi G ^(H), its projection /x onto lies in (/(II^); 

2. Any Gihhs measure fi G ^(H) is uniquely determined hy its projection onto 

3. Given any fi G ^(II^), there exists fi G Q(J\) such that the projection of ji 
onto is equal to fi. 

Proof. Let fi G ^(n) and let fi be its projection onto . Then, for any 
[/ G V and r G [recall the definition of given in (3.2)], 

^(au = 'ru\(7v\u = Tv\u) = K^v = Tv\^v\u = tv\u) 

(3.9) 

= H K^ = ^vv\^v\u = 'rv\u)- 

Given any rj G , fi{a = tv'>]\^v\U = Tv\u) is equal to 

f^^^UU[EndU] ^ '^^^Endu\^[V\U]u[E\dU] ^'^v\uVE\du) 
^ f^i^E\dU=^E\du\^V\U = 'Tv\u)- 

Note that U VJ[E r\ dU] G V, since C/ G V and the tree is a locally finite 
graph. Moreover, V\[U\J[Er\ dU]] = [V\U]U[E\ dU]. Therefore, using 
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the relation d\U VJ[Er\ dU]] = dU '^V\U proved in Lemma 3.5, along with 
the fact that fi satisfies (3.5), we see that 

f^^^uu[Endu] ^ '^u^Endu\^[v\u]u[E\du] ^'^vxuVExau^ 

The last two displays, when substituted back into (3.9), show that for every 
U eV and r G i^c, 

^i{(^u = Tu\<yv\u = rv\u)= Yl ^uuiEndu]('^du;^uuiEnau]='^uVEndu) 

^ i^(^E\dU = ^E\au\^v\u = Tv\u) 
= X! '^uu[EndU](^du;(^uu[EndU]^'^u'n') 

^ H K^E\dU=^"\^V\U = Tv\u) 
r^"(zSE\dU 

= X! '^uvj[EndU]^'^du'-,<^uyj^Er\dU]^'^u^') 

= ^uu[EndU](^du;^u = Tu), 

where the second equality uses the fact that, for any [/ € V, 5^ = S^^^^ x 
gE\dU ^ and the third equality follows because of the fact that, since /x is a 
probability measure concentrated on S'^^^ , the second summation in the 
previous line is equal to 1. The last display, when combined with Lemma 
3.6, then shows that for every U £V and /i a.s. r G 5^, 

fi{au = TuWv\u = Tv\u) = '^<t>,u{Tdu](^u = tu), 

which, by Definition 3.4, implies that fi S Q{Il^), thus establishing the first 
property of the theorem. 

To prove the second property, let fi G ^(H) and let /x be its projection 
onto . Then, for ?7 G V and r G , fi{au = tu) is equal to 

X! K^AU = Tu'n)= K^^unE = TuE'nounE\^^unv = TUv'nounv^ 

X K^Aunv = TVvVaunv) 
= K^AunE = TuE'naunE\^Aur\V = TUvflaucw) 

X niaAunv = TUv'naurw)-' 
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where the last equahty uses the fact that /i is a projection of jl onto S . 
Combining the above display with the fact that fl satisfies (3.5), and using 
the relation (9(A[/ n -B) = At/ n V proved in Lemma 3.5, we see that 

fi[au = Tu)= ^ T^AunEiruvVQunv'-'^'^unE = TuE'n9unE^ 

X ^J'{crAUnV = nivVdunv)- 

Since the last relation holds for every U £V, given /i, it uniquely character- 
izes fi. 

For the last property, suppose // G ^(II^). Let 9 be the probability distri- 
bution on given by 

g(i) = for i = 0,l,...,CE. 

Then, define fi as follows: to generate a sample a from /i, first pick a con- 
figuration a according to ^i, and then, independently on each edge xy of E, 
choose a configuration axy on Sc^ according to the distribution 9 condi- 
tioned on the constraint that dxy + cr^ + (Jy < C . As a direct consequence 
of this definition, we see that ft is supported on and the projection of 
fl onto is equal to /i. Using arguments similar to those given above to 
prove properties 1 and 2, it is easy to verify that, in fact, fi G ^(H). We leave 
the details to the reader. □ 

Corollary 3.7. Let II and be as defined in Theorem 3.1. Then, 
there is a one-to-one correspondence between the sets ^(H) and QiTlff^). 

Proof. Let x be the mapping that takes fi G G{tl) to its projection 
onto . Then, property 1 of Theorem 3.1 shows that x maps ^(n) to 
^(n,^), property 2 shows that x is one-to-one and property 3 shows that the 
mapping is onto. Thus, x defines a one-to-one correspondence between the 
two sets of Gibbs measures. □ 

As a consequence of Corollary 3.7, the analysis of phase transitions for 
Gibbs measures that correspond to the simple product-form specification 11 
that is Markov with respect to the more complicated graph G is reduced to 
the analysis of Gibbs measures corresponding to the more complicated spec- 
ification n,^, which is Markov with respect to the simpler tree graph G. The 
advantage of this reduction is that Gibbs measures that are Markov with 
respect to trees are usually much easier to analyze than Gibbs measures that 
are Markov with respect to more complicated graphs. However, it is worth 
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emphasizing that, while the original Markov specification fl is of product- 
form, the associated Markov specification 11^ is no longer of product-form. 
In addition, the set of allowable configurations in cannot be represented 
as the set of homomorphisms from the tree to a constraint graph, as was 
done in [5]. Hence, this problem does not fall directly within the frame- 
work considered in [5]. Instead, we use the full generality of the results in 
[28, 29], which consider Gibbs measures associated with general Markov 
specifications on trees, to prove Theorem 1.1, which shows that the analysis 
of Gibbs measures for the controlled unicast-multicast model can be further 
reduced to the study of recursions of the random field map given in 
(1.6). These recursions are analyzed for the case Cv = 1 in Section 4. 

Proof of Theorem 1.1. Let be the set of vectors -0 G M^^"*"^ 
having first component i/'q = 1. Let the mapping F'^ : — > "i^^ be given by 

Define i/''-'^^ € R^^"*"^ to be the vector whose first component is equal to 1 
and remaining components are equal to zero, that is, V^*^^ = (IiOj-'-jO); 
and let = F''^(t/>(")) be the iterates of under the map F<^. Then, 

Theorem 4.1 of [29] states that the Markov specification 11,^ has a unique 
Gibbs measure if and only if the sequence of iterates tp^^\n>0, converges 
to a limit and, in this case, the limit tp* must be the unique fixed point of 
the map F'^'. 

Consider the map A:"^^^ R^^' defined by 

A^W = ^fc^^'^'Vfc for A: = l,...,Cy, 

and let A~ :]R_|_^ — ^'(^ be the map given by 

r 1, for A; = 0, 

^^^ = W^''-''k., forfc = l,...,Cv, 

where is used to denote the fcth component of A{(j)). Then, clearly, 

AoA~^ and A~^ oA are the identity mappings on M.^^ and respectively. 
We now claim that, for any -0 E "^^^ 

(3.10) F*(V') = A~^<l>'^'^{AiP). 

It is easy to see that the theorem follows from this claim. Indeed, let be 
the zero vector in and note that = A{tp^'^^). Moreover, the relation 
(3.10) shows that ■0* is a fixed point of F"^ if and only if A{'ip*) is a fixed 
point of <^'^'^ and, by induction, that A[F'>']"- (ip) = [<^'^'^]"-{Aip). From this. 
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it immediately follows that ^ ^* if and only if [$'^'^]"(0) 

j4('0*). With this equivalence, the theorem is a direct consequence of the 
result from Theorem 4.1 of [29] quoted above. 

So it only remains to prove the claim (3.10). Suppose ip S Then, 
substituting the definition of 4> from (3.8), we have, for i = 1, . . . , Cy, 



,min( 



^mm{C-i,Cv) / \l/(q+l) x^'^in{C-{i+j),CE) x „- 



^2^0 

Y^min(C,CE) > Y^min(C-fc,Cv) 1/(9+1) ; 

^mm{C-i,CE) X ^min(C-(j+fc),Cv) , l/(g+l) 

,9/(9+1) f Z^fc=o ^fc Z^j=o '^i 

Y-Cb >, ^mm(C-fc,Cv) ,1/(9+1) 
-l/{q+l] (l^k=0 M[i- + l^j = l ^j{iP)\\ 
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^ E^foA.[i + E7=if''-''''^^^,(V')] ) 

where the third equality is obtained by exchanging the order of summation, 
the fourth equality uses the fact that C^; < C (which holds by assumption) 
to replace min(C£;,C) by Ce and the fifth equality holds because, for every 
/c = 0, . . . , C — min(C — i — k, Cy) > and, likewise, for every = 0, . . . , C, 
min(C — k,Cv) > 0. This establishes (3.10) and, therefore, completes the 
proof of the theorem. □ 



4. Analysis for the case Cy = 1? Ce = C- In Section 3 we established a 
relation between Gibbs measures for the controlled unicast-multicast model 
on the {q + l)-regular tree and recursions of the random field map 
In this section we analyze this recursion for the case Cy = 1 and Ce equal 
to some integer C > 2. (Note that the restriction C > 2 is without loss of 
generality, since, as mentioned earlier, when Ce = Cy = 1, the model is 
equivalent to a pure multicast model with a modified weight, which is well 
understood [18, 27, 28, 29].) When Cy = 1, u = can be parametrized 

by just the scalar u, and ^'^''^ = <I>'^''^:R+ reduces to the following 
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one-dimensional map: 



e£o^ Ta=o \ \ 



Ac + ^Ac~i 

where one should recall that A„ = J27=o "^i- Throughout this section, we will 
assume that the weight vector A satisfies condition (1.7) stated in Assump- 
tion 1.1. In Lemma 4.1, we now verify this condition for the weight vectors 
defined in (1.4) and (1.5). 

Lemma 4.1. Given any positive scalar \, Assumption 1.1 is satisfied by 
the weight vectors A defined in terms of X by (1.4) and (1.5). 

Proof. Fix C > 2. We first consider the case of the "Poisson" weight 
vector defined by (1.4), with Aj = X^/il. In this case, 

A^_i - AcAc-2 = Ac_i yAc^2 + ^c-l)l ) ~ ^ a)^^'^ 



(C-1)! C! 
A^^-i ez,V A^-1 A^A^ 

+ E 



(C-1)! ^^V(C-l)!(i + l)! C\i\ 

Now, the coefficient for A*^"*"* is positive for 0<i<C — 2if and only if 

1 1 11 

>0, 



(C-l)!(i + l)! C\i\ 

and this reduces to the condition that C > i -|- 1 , which certainly holds for 
< i < C — 2. Thus, each term is greater than zero, and (1.7) follows. 

Now, suppose the weight vector is defined by (1.5), so that Aj = A*. Then, 
simple algebraic manipulations show that 

(4.1) Al_,-AcAc-2 = X^-\ 

which is positive, since A > 0. Thus, once again, (1.7) is established. □ 



In Section 4.1 we fully characterize the phase transition surface (namely 
the boundary in parameter space between the regions of uniqueness and 
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multiplicity of Gibbs measures) for this model and show that the phase 
transition, whenever it occurs, is nonmonotone in the parameter v (for fixed 
g, C and A). In Section 4.2 we study monotonicity of the phase transitions 
with respect to the other parameters q and A. 

4.1. Characterization of the phase transition surface. In Lemma 4.2 we 
show that the map ^"^'^ always has a unique fixed point. 

Lemma 4.2. Suppose A satisfies Assumption 1.1 and i' > 0. Then the 
following two properties hold: 

1. The map ^'^''^ has exactly one fixed point, ^* = ^*{u,X), in the interval 
[0,(X)). 

2. For a fixed A, the range ^ £ (0, oo)} is equal to (0,00). More- 
over, there is a one-to-one correspondence between v and the correspond- 
ing fixed point given hy v = J{£,*), where J is the mapping defined in 
(1.8). 



Proof. Observe that d^'^' {^)/d^ is equal to 

'Ac^i+^Ac^2Y AcAc-2-Ac-i 



C-l 



uq 



Ac + eAc-i / (Ac_i + eAc_2)(Ac + eAc- 
which, in turn, shows that 



(4.2) ^ = $'^.^(£). 

^ ^ ^^^(Ac_i + eAc~2)(Ac + eAc;_i) 

Due to (1.7) and the fact that A^, <I>'^'''*(^) > 0, the last expression is negative. 
Furthermore, $^'^(0) = u{Ac-i/Ac)'^ > and lim^^oo ^'''^(O = ^{^^0-2/ 
Ac-iY- Hence, since <&'^'^(^) is continuous and strictly decreasing in ^, with 
^^'^{0 > for all ^ > 0, (^"'^ has a unique fixed point. 

To show the second property, fix A and note that v satisfies the implicit 
equation 

where ^* = is the unique fixed point of <^^'^. Evaluating the derivative 
of J at the fixed point, we obtain 



NONMONOTONICITY OF PHASE TRANSITIONS 



23 



Ac-i+eAc-2/ (Ac-i + eAc-2)2 

" e '^'^ (Ac + eAc-i)(Ac-i + eAc;-2) ■ 

By Assumption 1.1, A^_]^ — AcA(7_2 > 0, and so this derivative is positive. 
Thus, by the imphcit function theorem, there is a one-to-one correspondence 
between u and Furthermore, it is easy to see from the map that 

^ as ^ and — > oo as z/ ^ oo, so that the range of € 
(0,cx)) is (0,oo). □ 

In order to investigate when iterates of the map ^'^'■'^ starting at converge 
to the unique fixed point , we wih need the following dynamical systems 
result. Given an interval / C M, a fixed point x* of a map F :! ^ I is said 
to be globally stable if, for all x £ I, F'^{x) — > x* (where represents the 
nth iteration of the map F) . 

Lemma 4.3. Let F:I^I be a thrice continuously differentiable map 
on a bounded interval / C M that has a unique fixed point x* . If F has a 
negative Schwarzian derivative, that is, if SF{x) < for all x € I, where 



(fF/dx^ 3f(PF/d. 
dF/dx ~2\dF/d: 

then x* is globally stable if and only if 



dx ' 



< 1. 



Moreover, if F is decreasing on I, then global stability of the fixed point x* 
is equivalent to the condition that F^''{1) x* , where I is the left endpoint of 
the interval I . 



Proof. The proof of the first statement is a consequence of standard 
results concerning one-dimensional dynamics (see, e.g., page 158 of [9]). The 
fact that F is monotone decreasing implies that FoF is monotone increasing 
and that F{x) < F{1) for every x G /, so that we can assume, without loss 
of generality, that F(l) = r is the right endpoint of the interval /. Thus, for 
any x £ I, we have 

f2'^(/) < < F^^'iFil)) = F^^'+^l) 

and 
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which, together, show that F^{x) x* , for all x € / if and only if F^[l) 
X*, as desired. □ 

We now show that the map ^'^''^ satisfies the conditions of Lemma 4.3. 

Lemma 4.4. Suppose Assumption 1.1 is satisfied and ^* = ^*{u,X) is 
the unique fixed point of the map Then, [<I>'^'''*]"(0) — > ^* if and only if 

Proof. First, note that, by definition of the random field map, it fol- 
lows that, for any G [0,oo), <I>''''^(,^) G [0, i^] and, therefore, we can assume, 
without loss of generality, that the map ^'^'^ is defined on a bounded interval 
/. As a consequence of Lemma 4.3, it is sufficient to show that ^"^'"^ is thrice 
continuously differentiable, decreasing and has a negative Schwarzian deriva- 
tive (with respect to ^). It is easy to see that ^'^'^ is thrice differentiable and 
in fact satisfies 



9? ^"'(Ac-i+eAc-2)(Ac + ^Ac-i)' 

g2^.,A(^^ ^ d^'^'^iO jq - l)AcAc_2 - (g + l)A^_i - 2^Ac_iAc_2 
de ~ (Ac-_i + eAc-2)(Ac + eAc-i) 

de ~ (Ac-i + eAc-2)2(Ac; + CAc-i)2' 

where G(0 = [{q - - (q + l)A^„i - 2CAc-iAc-2][iq - 2) Ac x 

Ac7-2 - {q + '2)^c-i - 4eAc-iAc-2] - 2Ac-iAc-2(Ac-i + eAc-2)(Ac + 
^Ac„i). The fact that ^'^'"^ is monotone decreasing then follows from the 
first equality above, along with Assumption 1.1. Tedious algebraic manip- 
ulations of the quantities above, also show that the Schwarzian derivative 
being negative reduces to the condition that 

0<^iq^-l){AcAc-2-Al-if, 
which holds trivially. □ 

We now derive some elementary consequences of Condition A introduced 
in Definition 1.1 that will be useful for the proof of Theorem 4.1. 



Lemma 4.5. Suppose Assumption 1.1 holds. Then, the following state- 
ments are true: 
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1. // Condition A holds for some {q,C,X), then 

(4.3) (1 + g)AcAc-2 + (1 - '?)A^-i < 0. 

2. Condition A holds for some {q, C, A) if and only if the quadratic Q defined 
in (1.9) has two distinct real roots < a- < < oo, with Q{a) < for 
a G (a_, Q;+). 

3. The relation in Condition A holds with equality if and only if Q has a 
positive real double root a- = a+. 

Proof. For the first property, note that if Condition A holds, then 

qiq + l)AcAc-2 <{q + lfAcAc-2 < (q - l)^A^„i < q{q - l)A^_i , 

so that {q + l)AcAc^2 < {q - l)Ac'_i, namely (4.3), holds. 

For the second property, recall that the quadratic Q has two distinct real 
roots if and only if its discriminant is strictly positive: 

[(1 - q)Al_^ + (1 + q)AcAc^2? - 4AcA^„iAc_2 > 0. 

Thus, to establish the second property, it suffices to show that the last 
display is equivalent to Condition A. Indeed, using the elementary equality 
(a + 6)^ — (a — 6)^ = Aab, the last display can be seen to be equivalent to 

[Ac7Ac-2 + A^_i + q{AcAc-2 - ^c~i)f 

> {AcAc^2 + ^c-i? - (AcAc7-2 - Ac-l)^ 

which reduces, after straightforward algebraic manipulation, to the condition 
that 

{q^ + l)(AcAc7-2 - A^_i)2 + 2q{AcAc-2 + Al_^){AcAc-2 - ^c^i) > 0- 

By Assumption 1.1, AcAc_2 — A^_]^ < 0, so that this reduces further to the 
condition that 

{q^ + l)(AcAc-2 - Ac~i) + 2g(AcAc7-2 + Ac-i) < 0, 
which can be rewritten as Condition A: 

iq+lfAcAc-2<{q-lfAl-i. 
The minimum of Q is attained at the point 

-[(1 + q)AcAc-2 + (1 - g)A^„i]/[2Ac™iAc-2], 

which is positive, by (4.3). Since Q{0) > 0, any real roots that exist must be 
strictly positive, with Q{a) < for q € (a_,a+). 

For the third property, the arguments above also show that, if the relation 
in Condition A holds with equality, then the discriminant is zero and the 
double root of the quadratic Q occurs at {q — l)Ac_i/Ac„2 > 0. □ 
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We are now in a position to give a precise characterization of the phase 
transition region. Recall the set A and Condition A introduced in Definition 
1.1 and the function J defined in (1.8). Moreover, as in Lemma 4.5, let a_ 
and a+ be the distinct positive roots of the quadratic Q and let v- = J(a_) 
and z/_|„ = J(a+). 

Theorem 4.1. Suppose Assumption 1.1 holds. Let ^* =^*(z^. A) be the 
unique fixed point of the map Then, \d^{(,*)/d^\ < I if and only if 

either {q, C, A) ^ A, or (q, C,X) G A and i/ G (0, U [z/+, oo). 

Proof. Substituting the fact that the fixed point satisfies ^* = <I>'^'-^(^*) 
into (4.2), we see that 



Cq{AcAc-2-A-l_,) 



(Ac-i + eAc-2)(Ac + eAc-i) 
e(7(Ac-i - AcAc-2) 



(Ac7_i + eAc-2)(Ac + eAc-i) ' 

where (1.7) was used to obtain the last equality. Therefore, \d^{^*)/d^\ < 1 
if and only if 

g*g(A^_i - AcAc_2) ^ ^ 

(Ac-i+eAc;-2)(Ac + eAc-i) " ' 
from which we conclude that 

(4.4) \d^{e)/dc\<i ^ Q(r)>o, 

where Q is the quadratic given by (1.9). 

First, suppose that ((/,C, A) ^ A. Then, by Lemma 4.5, it follows that the 
quadratic function Q has no distinct real roots. Since Q{0) = AcAc-i > 0, 
this implies that Q{a) > 0, for all a G R+ and so, in particular, Q{d'^, ^)) ^ 
0, for all v G (0,oo). Now suppose that (g,C, A) G A, so that Condition A 
is satisfied. Then, by the second property of Lemma 4.5, it follows that 
Q has two distinct positive real roots a_ and a+, and Q{a) is negative 
precisely when a G (a_,a+). By the second property of Lemma 4.2, every 

a G (0,a_] U [a+,oo) corresponds to a fixed point of the map with 
u = J{a) G (0, J(a_)] U [J(q+), oo) and, thus, (5(^*(j^, A)) > for precisely 
these values of v. The theorem then follows from the above arguments and 
(4.4). □ 

The proof of Theorem 1.2 now follows. 



Proof of Theorem 1.2. The proof follows as a direct consequence of 
Lemma 4.4, Theorems 1.1 and 4.1. □ 
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4.2. Monotonicity of phase transitions with respect to q,v, A. Theorem 
1.2 provides an exphcit characterization of the region in the ((?,C, A,z^) pa- 
rameter space where there are multiple Gibbs measures. Note that, in par- 
ticular, this result shows that if there is a phase transition at all for some 
(g, C, A, v), then this phase transition is nonmonotone in z^, in the sense that, 
for the same ((7,C, A) and every sufficiently small and sufficiently large z^, 
there exists a unique Gibbs measure. The numerical examples in Section 5 
show that phase transitions do indeed occur — in the sense that for a large 
class of weights, the set A\s a, nontrivial subset of the parameter space — and 
also illustrate the nonmonotonicity with respect to v. Some intuition behind 
why this nonmonotonicity occurs is provided in Section 6. 

We now establish some additional monotonicity properties of the phase 
transition. The first result. Theorem 4.2, shows that the phase transition 
surface is, in a sense, monotone in q, the degree of the tree. 

Theorem 4.2. Suppose Assumption 1.1 is satisfied. If {q,C,X) (z A, 
then {q + 1,C, A) € A. 

Proof. First, note that Condition A holds for q + 1 if and only if {q + 
2)^AcAc_2 < q^Is^_i, which, in turn, holds if and only if 

{q + lfKc^c-2 + {2q + 3)AcAc-2 

(4.5) 

<{q-lfKl_^ + {2q-l)Al_^. 

Now, suppose Condition A holds for q. This implies that {q + l)^AcAc_2 < 
(g — 1)^A^_-^ and, due to the first property of Lemma 4.5, that {2q + 2)Ac x 
Ac_2 < (2g — 2)A^_-^. In addition, due to Assumption (1.1), we also have 
AcAc_2 < A^_^. When combined, this shows that (4.5) holds and, therefore, 

that Condition A also holds for g + 1, or, equivalently, that {q + 1, C, A) G A. 

□ 

We now show that, for the weights in (1.4), the set A is also monotone 
with respect to the parameter A. 

Theorem 4.3. Suppose that X is a weight vector defined as in (1.4) for 
some A > 0. Then, given any weight vector X' defined as in (1.4) but with X 
replaced by X' > X, {q, C,X) eA^ {q, C, A') G A. 

Proof. We first show that {q,C,X) e A if and only if F{q,C,X) > 0, 
where 



F{q, C, A) = (1 + qf[Xc-iAc-i - XcAc-2] - 4<?A^- 
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Indeed, by Definition 1.1, {q,C, X) G A and only if 

(l + (?)2Ac7Ac-2<(g-l)'A^_i 

^ {q + l)2(Ac-i + Ac)(Ac7-i - Ac-i) <iq- l)'Ac-i 

^ (1 + qf[Xc~iAc-i - Xc^c-2] - 4gA^_i > 0. 

Therefore, in order to prove the theorem, it suffices to show that if F{q, C, A) > 0, 
then dF{q,C, X)/dX > 0, since the latter implies that, for any A' > A, F{q,C, 
A') > 0, which, in turn, means that (g, C, A') € A. 

Let H{i) be the coefficients in the expansion of F{q, C, A) as a polynomial 
in A: F{q, C, A) = X]?=(r^ H{i)X^ . We now claim that the following assertion 
is true: 

Claim. There exists C -1 <i* <2C -2 such that H(i) < for all i < 
and H{i) > for all i>i*. 

We defer the proof of the claim and, instead, first show that the the- 
orem follows from this claim. Indeed, if the claim holds, then for i > i* , 
iH{i)X'/i* > H{i)X' > and, for i < i*, > iH{i)X''/i* > H{i)X\ Therefore, 
we have 

XdjW^ = A .i^(^)A-l = V + -H^' 

i* dX i* ^ ^ ^ ^ i* ^ ^ i* ^ ' 

1=1 1=1 1=1* 

2C~2 

> J2 H{i)X' = F{q,C,X)-H{0) 
1=1 

>F{q,C,X)>0, 

where the second-to-last inequality follows because H(0) < 0. Thus, when- 
ever F{q, C, A) > 0, we also have dF{q, C, X)/dX> 0, which (as argued above) 
implies that A is monotone in A. 

Thus, it only remains to prove the claim. Note first that, for < i < C — 1, 
H(i) < 0, since A* appears only in the —AqAQ_-^ term. Now, if F{q, C, A) > 0, 
then it cannot be that H(i) < for all < i < 2C — 2 and, therefore, there 
must exist a minimal C — 1 <i* < 2C — 2 such that H{i*) > and, for all 
i <i* , H[i) < 0. Thus, to prove the claim, it suffices to show that H{i) > 
for all i>i*. 

Now, for C — 1 < i < 2C — 2, observe that the coefficient H{i) for A* is 
given by 



H{i) = {l + qf 



1 ^ . C-l 



{C -l)\{i-C + C\{i-C)\ 



k=i-C+l 



kl {i-k)\ 
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(1 + qf- ((C!(z -C + 1)!)/(2C - I - l))4gEEr/,c+i(lA!)l/(^ - k)\ 
{C\{i-C + l)\)/{2C -i-l) 

It is easy to see that then H{i) > if and only if 

-~.r.r. ^ 4C!(i-C + l)! t,' 1 1 , 

Now, note that if 

(4.6) G{i + l)<G{i) for i = C-l,...,2C-3, 

then H{i*) > 0, which imphes q > G{i*), which, in turn, imphes that, for 
all i > i*, q> G{i) and, therefore, H[i) > 0. Thus, in order to establish the 
claim, it suffices to prove (4.6). 

To establish (4.6), we will find it convenient to first change variables, 
namely, to set j = 2G — 1 — i and m = C — k for A; = 1, . . . , C — 1. Then, for 
i = l,...,C, 



4C!(C-i)! ^ _J- 1 

3 ^,(C-m)!(C-j-l + m)! 

4C^ (C-l)!(C7-j)! 



3 ^^{G-m)\{C-j-l + m)\ 



where we define 



4C 

— Er(j,m)-2, 

7 — 

m=l 



{C-\)\{C-i)\ 

T{j,m) - 



{C-my.{G-j-l + m)l' 

We now need to show that G{j) < G{j + 1) for j = 1, . . . , C — 2. 

We first summarize some properties of T[j,m) that will be used in the 
sequel. It is easy to see that T(j, m) = T{j, j — m + 1) and, after some manip- 
ulation, that T(j, m + 1) = {G — m)T(j,m)/{C — j + m). Thus, T(j, m + 1) > 
T(j, m) if and only if (G — m)/{C — j + m) > 1, that is, if and only if m < j'/2. 
Similarly, T{j,m + 1) < T{j,m) for m > j/2, and T{j,m + 1) = T{j,m) 
for m = j/2. Consequently, for a given value of j, T{j,m) has its maxi- 
mum value at m = j/2 when j is even, in which case it takes the value 
(C - 1)!(C - j)!/[(C - i/2)!(C - j/2 - 1)!], whereas,when j is odd, the 
maximum (C - 1)\{G - i)!/[(C - (j + 1)/2)!(C - (j + l)/2)!] is achieved 
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at m = (j + l)/2. It is also straightforward to check that T{j + l,m) = 
(C — J — 1 + m)T{j, m)/{C — j) > T{j, m), since m is positive. 

Now, let us consider G{j) again. If j is even, we split the sum as follows: 

■ i/2 j/2 

^0' ^(•?' J - m + 1) 

.m=l m=l 

where we have used the fact that T[j,m) = T{j,j — m + 1). Then, again 
making repeated use of the properties of T(j, m) listed above, note that 

AC 



J + 
AC 

7+1 



m=l 

ri/2 



^r(j + i,m) + r(j + i,j/2 + i) 

.m=l 

i/2 

+ Y,T{j + lJ-m + 2) 



m=l 



> 



AC 

J+T 



i/2 



.m=l 



AC 



Now, 



j + 1 

r(j + i,jy2 + i) 



5]r(j,m) + r(j + i,j72 + i) 
+ 5^ro-,j-m + i) 

m=l 

^r(j,m) + r(j + i,jy2 + i) 

(C-l)!(C-j-l)! 



.m=l 



{c-j/2-iy.{c-j/2-iy. 

C-j/2, 



-T{j, j/2) 



c-j 

>T{j,j/2). 

On the other hand, T{j,j/2) is the maximal value of T{j,m), and so 
T{j + l,j/2 + l)>T{j,m) 

j 

=^ jT{j + i,j/2+i)>j2nj, 

m=l 



r(i + i,j72 + i)>-^r(j,m). 

7 — 

•' m=l 
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Hence, 



G{j + l)> 



AC 

JTT 



J2T{j,m) + WT{j, 

.m=l •' m=l 



m] 



AC 

= —Y,T{j,m)-2 
1 — 

■' m=l 

= G{j). 

If j is odd, a similar argument gives the result. We omit the details. Thus, 
for j = 1, . . . , C — 1, we have shown that G{j + 1) > G{j) with equality if and 
only if j = 1. Returning to our original variable i = 2C — 1 — j, we therefore 
have G{i + 1) < G{i) for C < i < 2C - 2 and G{2C - 2) = G(2C - 3), which 
completes the proof of (4.6) and, therefore, of the claim and theorem. □ 



Remark. It is worthwhile to note that the set A is not monotone in A (in 
the same sense as specified in Theorem 4.3) when the weight vector A is given 
by (1.5) for some A > 0. To see this, consider, for instance, the case C = 2. 
Then, Condition A reduces to the inequality (l + g)^A — 40^(1 + A)^ > 0, which 
clearly fails to hold for A sufficiently large. In particular, when (7 = 14 (with 
C = 2), Condition A holds only for the weight vectors A associated with A G 
(49/56,64/56). In other words, for q = 14, C = 2, there are multiple Gibbs 
measures for the parameters (g, C, A, i^) if and only if A € (49/56, 64/56) and 
1/ lies in a nonempty open interval (z^_(A), z^_|_(A)), while, for all other values, 
there is a unique Gibbs measure. 

5. Numerical examples. In this section, we illustrate our results with 
some examples. 

An example of nonmonotonicity of the phase transition with respect to 
the multicast arrival rate v can be seen in Figure 1, which plots the blocking 
probability for multicast calls when Ce = C = 2, Cy = 1, A = 0.75 and q = 
10, with weights Aj = A*/i!. It can be clearly seen that there is an interval 
(j>_,i/+) within which the fixed point is unstable, while outside this interval 
it is stable. This should be contrasted with the case where C = 1 and Cy = 1, 
where it is known [26] that the phase transition is monotone in v. In this 
example, we see that, as a result of the control, a difference in blocking 
probabilities in the even and odd sublattices of the tree, though it exists, 
is not significant and, hence, the unfairness due to heterogeneous blocking 
probabilities is not severe. 

Figure 2 plots A against and v+ when Ce = C = 2, Cy = 1 and q = 6. 
Note that, for these parameter values, there is no phase transition when 
A < 6. The region enclosed within the curve is the parameter region where 
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the Gibbs measure is not unique. One very interesting feature that is not 
immediately apparent from this plot is that, for A G (6.000,6.103), de- 
creases as A increases. In other words, the point at which multiple Gibbs 
measures appear may decrease as A increases. This contrasts with the case 
C = 1, where an increasing A always increases the point at which multiple 
Gibbs measures appear. 

Figure 3 is a similar plot for Ce = C = 2 and Cy = 1, but this time with 
g = 14 and weights Aj = A*. This illustrates the result that, for these weights, 
the occurrence of multiple Gibbs measures is not monotone in A; multiple 
Gibbs measures only occur for A G (49/56,64/56). 

The phenomenon of nonmonotonicity appears not to be restricted to the 
controlled model, with Cy < C. We conclude this section with the plot in 
Figure 4 of the multicast blocking probabilities when Ce = Cy = C = 2, 
q= 10 and A = 0.72 for v € (10, 60), which illustrates this point. 



6. Concluding remarks. In this paper, we analyzed an idealized model 
of multicasting on a regular tree network with a simple control that limits 
the number of multicast calls that can be centered at any node. With this 
control in place, we showed that a phase transition may exist in the infinite 
network and that it is nonmonotone in z^, the arrival rate of the multicast 
calls. 



s 



i 



1 ^ 




Fig. 1. Multicast blocking probability for weights \i — /i\ (Ce = C = 2, q = 10, Cv = ^ 
X^O.75). 
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Fig. 3. The phase transition region, C = 2, q = 14, Cv = 1, Ai = A*. 
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Fig. 4. Multicast blocking probability for a model without controls, Cv = Ce = C = 2, 
g= 10, = A = 0.72. 

It is in general a nontrivial task to identify when a phase transition is 
monotone with respect to the weights. For example, the identification of 
graphs for which the phase transition of even the hard-core model is mono- 
tone remains an open problem and, in particular, this question has not been 
resolved even for the d-dimensional lattice (see [16] for some partial results in 
this direction). Nonmonotonicity of phase transitions for a hard core model 
on a tree network with added hanging links was observed by Brightwell, 
Haggstrom and Winkler [4]. In this work, we presented an example of non- 
monotonicity with respect to a parameter that has arisen naturally in the 
context of controls on multicast calls in a regular communications network. 

In order to gain insight into the nature of the phase transitions studied 
here, it is useful to examine the different ways in which the two types of 
calls (i.e., unicast and multicast calls) can be packed into the network. A 
nonrigorous, but intuitive, explanation of the nonmonotonicity of the phase 
transition with respect to the multicast arrival rates is as follows. For sim- 
plicity we concentrate on the case when Ce = C = 2 (with Cy = !)• Then, 
for all sufficiently low z^, the network behaves more or less as it would in the 
absence of multicast calls, in which case (since unicast calls at different links 
do not interact with each other) there is clearly a unique Gibbs measure. 
This unique measure can be thought to roughly correspond to a packing 
of two unicast calls throughout the network, which is clearly homogeneous. 
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However, as v increases, there is increasing demand for capacity by multicast 
calls, leading to competition between multicast calls and unicast calls. As a 
consequence, the unicast calls are now restricted to occupy only one unit of 
capacity throughout the network, resulting in the multicast calls effectively 
experiencing a tree network with capacity 1. In this case, as is well known 
[18], there are two phases, with one corresponding to a packing of one mul- 
ticast call on each node on the even subtree, and the other to a packing 
of each multicast call on the odd subtree. As the multicast arrival rate in- 
creases yet further (relative to the unicast rate), the unicast calls are unable 
to capture even one unit of capacity throughout the network and, thus, one 
obtains once again a single homogeneous Gibbs measure corresponding to a 
packing of a multicast call at each node of the network. Similar reasoning 
leads us to the following conjecture, which is supported by numerical results 
(see Figure 4 of this paper and Section 5 of [21]). 

Conjecture. For the uncontrolled unicast-multicast model with Ce = 
Cy = C = 2, there exist phase transitions that are nonmonotone in the mul- 
ticast arrival rate v. 

Recent years have witnessed an explosion of research papers devoted to 
the analysis of phase transitions in statistical physics-type models in various 
fields, with the majority arising from problems in combinatorics and the 
theory of computing (see, e.g., [3, 6] and references therein). This work, along 
with [13, 17, 24] and [26], demonstrates that Gibbs measures associated with 
loss networks also exhibit several interesting phenomena that are worthy of 
further study. 
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